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C^l I Abstract. We suggest a theoretical picture for distributions of plastic 

deformations experienced by a sliding Charge Density Wave in the course 
of the conversion from the normal current at the contact to the collective one 
in the bulk. Several mechanisms of phase slips via creation and proliferation 

^ ■ of dislocations are compared. The results are applied to space resolved X- 

Tjj- ! ray [1,2], multi-contact [3,4] and optical [5] studies. Numerical simulations 

O ' are combined with model independent relations. 
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1. INTRODUCTION. 



Plastic deformations and flows are currently the focus of studies of Electronic 
Crystals (charge/spin density waves(DW) and interface 2d Wigner crystals) and 
^ • of vortex lattices. Thus, dislocations are supposed to participate in depinning, in 

^ I phase slips, in narrow band noise generation, and also arise in contact structures. 

P ■ Experimental studies show that the sliding state of DWs is essentially inhomo- 

geneous [1-6] with different effects near contacts and in the bulk. We present 
general scheme to describe inhomogeneous distributions of deformations (which 
^ , are the gradients (f' = q along the chain axe x of the phase (p), electric fields and 

^ I currents in the course of the current conversion. For a long enough sample we pre- 

dict a strong variation for the shift q{x) of the CDW wave number near contacts 
which flattens (nearly exponentially) towards the bulk, 5x — »• cxD.The bulk strain 
like q = q'f^i^x remains only when the conversion is not able to bring the partial 
currents to equilibrium with strains allowed at this expense. It happens either 
for size dominated configurations like in [3], or when the pinning of dislocations 
(DLs) hinders the conversion process as demonstrated by X-ray experiments [1]. 
We present numerical modeling for time dependent and stationary cases. 

2. MODEL INDEPENDENT RELATIONS. 

Generally in the CDW state there are two types a = i, e of normal carriers 
which share the metallic state DOS Np ( per chain area s). For intrinsic ones, 
a = i, the partition Np = f3iNp is closed below Tp, so that they are only excited 
across the Peierls gap. The extrinsic carriers are present in semimetalic DWs: 
NbSe^ and some SDWs. Their DOS partition PeNp remains open, /3i + /5e = 1- 
The potentials V^ experienced by the carriers and their chemical potentials /!„ 



are [7,8] 

F, = $ , \/, = $ + q/iTiN'p) ; /i, ^ \4 + 5nJ{p^N'^) 

where $ is the electric potential and 5na are the local imbalances between electron 
and hole concentrations for each type of carriers. The normalized normal and 
condensate densities, ptiT) and pdT) = 1 — Pi(T), change from pi = 1, pc = in 
the metallic state to pi = 0, pc = I at low T while Pe ~ 1 stays constant. Fast 
equilibration of quasi-particles leads to pi = pe = fJ^n and the normal current 
is jn = —{cFns/e^)dpn/dx whcre a^ = Uj + Ue is the normal conductivity. (All 
quantities are given per electron per chain, hence the conductivities acquire the 
factor s/e^.) 

We separate the total densities of the charge (5ntot and of the current jtot into 
their normal and condensate counterparts, 5nx,ot = ^^n + ^^c and jtot = jn + jc, 
with (5n„ = Sn^ + 6ni, 6nc = f'/'n' and jc = —(p/TT. (The renormalization 5nc — ^ 
Pc'-p' JT^ appears inherently [8] as an action of the potential l^.) The stress 2C/, 
per DW period, includes [7] three contributions: elastic, electronic and electric. 

U = (g/vr + 6ni)/N'p + $ , dU/dx = -T = Tire (1) 

The electronic term ~ 5ni in U is conjugated to the term ~ g in Vj. Both come 
from corresponding variations of the breathing energy 6niq/{nNp). As a variation 
of the energy over 6nc, the potential U plays a role of the condensate potential 

The second equation in (|^) states that the driving force of the DW sliding JF 
is equilibrated by the friction force —J^frc (jc) which coincides with electric field 
£ in equilibrium (g = 0, 6ni = 0). Thus the function J-'{jc) can be obtained 
from the sliding V — I characteristics for distant nonperturbative contacts. We 
always suppose that JF and jtot are above the sliding threshold values J-'t and 
jt = {.o'ns/e^)Tf The local electro-neutrality condition [9] implies that drin images 
q or vice versa: — (5n„ = drti. = q/ii thus linking X-ray [1,2] and optical [5] 
experiments. With the above relations we can express q and £ as: 

qg/T^N'p = Pn-U = r] £ = ^(jj - q'(3,g/T^NU g-' = (3, + pj p, (2) 

Here g = g{T) is the normalized elastic modulus with (3^ and pi characterizing 
the Coulomb screening of the CDW deformations by carriers. (In other theories 
only g = Pc has been exploited which is valid only near T^.) 

We see that only via extrinsic carriers the elastic stress enters the balance 
between the forces. Otherwise £ follows J^{jc) identically. The first relation in 
(^ shows that the CDW is deformed whenever the stress U does not coincide 
with the chemical potential of the normal carriers /i„. The quantity rj measures 
this mismatch, and hence characterizes an excess or a lack of normal carriers 
which is imaged directly by q. The exchange between the normal carriers and the 
condensate via the phase slips finally equilibrates pn and U, a process taking place 
via nucleation and growth of DLs. The asymptotic condition for the bulk is that 
all partial currents are stationary which implies the absence of current conversion, 
all types of carriers being in equilibrium, with equal chemical potentials hence 
?7 — >■ (3. Only hindering of phase slips can prevent leveling of rj, and hence q, to 
zero far from a contact. 

The population of normal carriers is controlled by the injection/extraction 

rates Uinj = ±jtotS{x — Xcnt) from the contacts and by the conversion rate 2TZ 



to/from the condensate: dun/dt = hn + j'n = i^inj — 27?.. Here TZ = lZ{rj,jc) is the 
rate of phase shps (per unit length per chain) which transform pairs of electrons 
to/from the DW periods. (See [10,11] for a review of microscopic mechanisms.) 
With the above relations we arrive at the two equations for two variables r] and 



wN'p/g + j^ - 27?(r7, jj + z/i^j = ^(j J - (eVsa„) {jt^t - 3c) = V (3) 

They provide a general scheme for the plastic sliding with time and space depen- 
dent current conversion processes. There are common elements with the modeling 
scheme used in [3]. The major difference is that our equations provide the flat 
asymptotics g ^ and also Eq. (2) suggests a different relation between the 
measured £ and the extracted jc- Next, all dynamics is generated by a nonlin- 
ear function J-'{jc) while in [3] its linear part (dominating at jc ^ jt only) was 
extracted, the rest being treated as an a posteriori correction. Further on, our 
form of g allows to study the semiconducting regime. 

The most sensitive element of any construction is the function TZ{r],jc) which 
may differ in various circumstances. This function is of great interest itself so 
that the above Eqs. can be utilized to recover it from measurements of q or S. 
Otherwise we can guess its form in plausible pictures of the current conversion 
to get q, S and to compare those with experiment. The HomoN scenario refers 
to an ideal host crystal, both in the bulk and at the surface, where only ho- 
mogeneous nucleation is present as a spontaneous thermal [12] or even quantum 
[13] supercritical fluctuation, so that T^hom oc exp[— r^o/l^l], valid at \ri\ <^ rjQ. In 
general the HomoN scenario is very unlikely in any physical system as requiring 
an extreme purity; also in our problem it cannot provide an observed fast contact 
variation of q. 

The HeteroN scenario refers to samples with a sufficiently large density of 
defects acting as nucleation centers for supercritical DLs, the simplest form for 
this case being 7?.hct ~ RpNpf]. Since Nptj/g is an over- saturation of the electronic 
concentration then gRp = t~^^ is nothing but the inverse life time of an excess 
normal carrier with respect to its absorption by a DL. In DWs the heterogeneous 
nucleation will be provided by any region where q varies over the cross-section 
which must be encircled by a sequence of DLs, one per a length of 2ti increment. 

Both scenarios are of the passive type, when the DW motion itself plays no 
role. A plausible active scenario emerges for a fast enough DW motion when the 
DLs are created by the DW sliding through bulk or surface defects: IZ = TZaiv^ Jc) 
such that 7?.a(0,jc) = "^0(^,0) = with IZa ~ RaNpTjjc as the simplest version. 
The origination of DLs at the crystallographic steps interrupting the DW sliding 
in a fraction of the sample cross-section is a most apparent manifestation. The 
origination of DLs pairs at defects in the bulk have been described in [10,14] 
within the theory of the nonlinear J-'frc{jc)- An extension of this theory to describe 
also TZ would bring the whole complex of sliding effects to one frame. 

At lowest |?7| the conversion bottleneck is a pinning of DLs by defects similar 
to vortex lines. It results either in a complete freezing of their motion or in a 
drastic reduction of their climb speed. Introducing r]t as such a threshold or a 
crossover we have the following guesses: i. collective pinning of DLs: TZ{i]) = at 
\ri\ < rjt; ii. local pinning of DLs: 7^(77) = RpNpf] exp{—rit/ \i]\) The last law for 
the DLs climb provides a logarithmic time decay of q observed near contacts in 
pulsed experiments of [1], similar to the Kim- Anderson law in superconductors. 



3. APPROXIMATE SOLUTIONS AND MODELING. 



For an illustration consider a linear sliding I — V: a^ = d{J-'/djc = cnst. 
(Actually it is true only asymptotically, for currents much above 2Jt as in ex- 
periments of [1] and even above 5Jt as in [3].) Then for the HeteroN model we 
arrive at the single linear Eq.: 



fj/g - r]"rl)a* + 2if]Tlp = ±5{x - x^nt)itot/N\ 



r 



D 



ATxe^N'p/s (4) 



where a* = cr„(Tc/(cr„ + Cc) and rjj ~ lA is the screening radius (by intrinsic 
carriers) in the metallic phase. This equation defines the characteristic time 
r = {2g7lp)~^ and length A = roi^na* /TZpY^"^ scales. For the local "phase slip 



voltage" Vps = $(a;, t) — xjtot/c we find 



Vps{x,t) = q{x,t) 



TfNi 



(5) 



The Eq.(H) shows that Vps indeed follows the strain q but with an unexpected 
coefficient, which may change the sign as a function of T or jtot at Pe = (^nlipn + 
(Tc). Notice that i— and e— carriers contribute oppositely to the contact electric 
field, see e.g. [8]. 

The apparent analytical solution of Eq. (^ allows for easy time dependent 
simulations e.g. for alternating pulses of [3]. Fig.l shows the evolution of q{x) 
with time after the change of the current polarity. We conclude that the non- 
monotonic behavior found in [3] is universal with respect to the conversion mech- 
anism. 
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Figure 1: Time dependent solutions (parameters from [3]). 
For stationary distributions the same linear model gives for Xcnt = io 

1 jtot A sinh(x/A) 

q = . 

4(7 o'n rjj cosh(a/A) 

4 



(6) 



Its fitting, corrected for a small cutoff rjt, has been proved [1] to be very good. The 
Eq. d^) shows that for a long enough sample there will be no linear gradient in the 
bulk. The carrier conversion will be complete and the CDW deformation reduced 
to zero. Otherwise for close, a < X, contacts of experiments [3] we find a lin- 
ear variation q' = jtot/{4:ganr%) corresponding to the essentially non-equilibrium 
regime with a negligible current conversion. For the nonlinear homoN model also 
from the analytical solution we find near one contact: q oc (In |a; — a|)~^. This 
is an extremely slow decrease which results in an incomplete carrier conversion 
and in a large strain along any realistic sample length, which contradicts to later 
X-ray [1] and multicontact [4,5] experiments. A reasonable fitting with this form 
of TZ [2,3,12] was due to its utilization well beyond the validity range |?7| -C r/o 
when actually one should change to power laws. 

The Fig. 2 shows the distribution 2q{x) for j = 2.13 jt plotted against the 
double shift data of dc experiments [1]. The x-axis gives the distance from the 
left contact at a; = —2mm to the mid point a; = 0. Taking from the data the 
contact value of 2gcnt = 12.6 x 10~''6* and assuming ^f = 1, we determine the 
contact overstress r^cnt = 1-2 meV. From our fit we find the cut off value r/t = 
0.24 meV in a good agreement with the value rjt = 0.25 meV derived from the 
persistent strain qt = 1.10~^ A~^ derived from pc data [1]. A comparison of the 
passive and active conversion mechanisms is given by the ratio Rp/Rajc ~ 1-2 
determined from our fit, hence the two contributions are comparable. 
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Figure 2: 2q{x) versus experimental points from [1]. 

In Fig. 3 we present calculations for the configuration of [3] within the same 
model and with the same parameters as for [1]. Now the distance is short L = 
670/im and the current is high jtot = ^jt- We cannot detect any visible difference 
between the simulations with the cut off rjt and without it. The pronounced 
bulk gradients observed for these short distances between contacts disappear 
completely when we increase this spacing 2.5 times which is just the rescaling 
from multi-contact to X-ray geometries. Hence in experiments [3] the whole 
length is affected by the contacts, the partial currents are far from equilibrium, 
conversion goes on over the whole sample but is far from being completed. 

Clearly in X-ray [1,2] and in multi-contact [3] experiments the two completely 
different regimes have been observed. Both experiments [1,2] recover the quenched 
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Figure 3: q{x) versus data from [3]. 

gradients of the same origin, but not the phase slip stress from contacts as it has 
been supposed. 

Our fitting allows to estimate a mean concentration (per sample cross-section) 
of the lattice faults Nfn which provide the active conversion. Indeed j'^ = 
2VNfitjc where V is the probability to create a D-line at one fault over-passing 
each DW period. Our hypothesis R = RaNpfjjc means that V{ri) = rj/rj* where 
the constant rj* is determined by the condition V{ri*) ~ 1. Hence r]* can be 
estimated as the highest energy payed to create a D-loop which corresponds to 
just one layer from the surface. From the characteristic energy of interplane in- 
teractions we expect rj* ~ lOOK so that Nfu ~ 10~^A~^. We have obtained 
the distance Lf^ = l/A^at ~ 10 /im which is less but comparable to the initial 
decay scale {q/q')cnt ~ 60 /xm and 30 times below the total characteristic length 
~ 300 /xm. We see that the current conversion is provided by about 10^ — 10^ 
defects with a maximal (near the contact) characteristic efficiency of few pro- 
cents: r]cnt/v* ~ 0.1. We conclude that D-loops of opposite signs predominantly 
annihilate each other by exchanging locally the normal currents and only some 
fraction survives to proliferate across the sample. 

4. DISCUSSION AND CONCLUSIONS. 

In sliding state the collective current which can reach its nominal value pro- 
vided that the equilibration between the normal and the collective components 
is completed. In contradiction to earlier theoretical expectations and interpre- 
tations of experiments, we show that this state is mainly undeformed. Namely, 
elastic deformations - either from the distributed stress [2,10] or from the phase 
slip tension at the contacts [3] are incompatible with an equilibrium with respect 
to the conversion. 

We conclude (see also [1]) that for a long enough sample the large g- values 
near the contacts fall nearly exponentially towards zero in accordance with the 
HeteroN scenario. The detailed fit requires for the contribution of the active 
conversion generated by the sliding itself. The large gradients observed in earlier 
multi-contact studies [3] can now be re-interpreted as due to the size effect: the 
distance between contacts is so short that only a small part of the applied normal 



current is converted, and the currents jc and j„ stay far from their equihbrium 
values. In this hmit the strain q{x) is close to linear with only small contact in- 
crements. In contrast, earlier x-ray experiments [2], performed on longer samples 
(4.5 mm) could not yet resolve the contact region and what was observed as a 
bulk gradient coincides with remnant distortion q'^ from pinned DLs as revealed 

by [1]. 

Till now the DW abilities to slide and to be deformed were supposed to be 
related. Our results show that they rather exclude each other. 
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